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CHAPTER 1 


INTRODUCTION 

Reducing the cost of microwave circuits goes hand in hand 
with reducing their size. In this respect, microstrip line, 
stripline, and suspended stripline have been essential in 
eliminating bulky waveguides and rigid coaxial lines in a 
great majority of microwave systems. Only in a few and very- 
demanding applications, such as high power transmission , or 
low-loss filtering, are waveguides still being used [7] . For 
realization of filters and oscillators, we need a frequency 
determining element which may be compatible with the 
microstripline, stripline or suspended stripline in size and 
with the conventional waveguide and coax resonators in 
performance, for miniaturization of microwave circuits. With 
the advent of high dielectric constant (cd) and high-Q 
dielectric resonators (DR) it is possible to meet above 
requirements and to provide an alternate solution for the 
conventional microwave resonators. Due to their wide variety 
of applications in microwave circuits, investigation of 
their circuit properties had been recognized as an important 
problem in early seventies [14-17] . These resonators offer 
various advantages in microwave field due to combination of 
properties they posses. Due to high dielectric constant, 
size of the dielectric resonator is smaller than the size of 
an empty cavity resonator operating at the same frequency, 
approximately by a factor of !/✓ ea. They fill a gap between 
wave guide and planer transmission lines like 
microstripline, stripline, suspended stripline, or inverted 
microstripline by providing high-Q and temperature stability 
comparable to those of invar cavity resonators along with an 
integrability reaching to planer transmission line 
resonators. They offer versatility and are adoptable to 
various microwave structures and coupling configurations. 
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Numerous analysis and synthesis methods exist for studying 
the behavior of dielectric resonators. Dielectric waveguide 
model (DWM) , Mode matching, Integral equation etc. are a few 
of them [7,14-18]. Recently, finite methods like finite 
element method and finite-difference time domain method, are 
also getting a great popularity with the availability of 
powerful computers [43] . 

In planer transmission lines , the microstripline has 
got a special importance for microwave integrated circuits, 
because of its compatibility and ease of manufacturing. But 
for higher frequencies especially in millimeter wave 
range (30-300 GHz), it faces the problems like higher losses, 
tight fabrication tolerances, handling fragility etc. and 
more seriously, it becomes very dispersive medium of 
transmission. The supplement to these problems are suspended 
stripline and inverted microstripline structures. These 
structures has lower losses, resulting in a higher Q-f actor 
than microstripline because a greater portion of the field 
exists in the air than that in microstripline. In addition, 
a broader line width is possible for a prescribed 
characteristic impedance, and the TEM mode propagation is 
more pronounced, with less dispersion . Thus, numerical 
quasi-TEM analysis can be applied to yield sufficiently 
accurate results up to the millimeter wave frequency 
range [8]. Hence, these suspended substrate structures are 
very useful for millimeter wave applications. 

Many references are available for getting the analysis 
of DR and the design parameters like resonant 
frequency [14-25] , coupling coefficient between two 
DRs[29,30], and coupling between DR and 
microstripline [31, 32, 42] in MIC environment, that can be 
used for low microwave frequencies. But, not much 
information are reported about these parameters of a DR in 
suspended- substrate environment. A need for improved 
technique for providing design solution of DR in this 
environment has generated the keen interest to find out 
efficient and accurate design procedures . 
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This thesis studies the circuit properties of DR in 
suspended substrate environment, specially for millimeter 
wave applications. Emphasis is given on cylindrical DR 
operating in TEOld mode, due to it's extensive use in 
microwave circuits and ease of exciting this mode. The main 
objective here is to develop suitable design criteria, so 
that dielectric resonators can be used as a circuit element 
in suspended substrate or MIC environment, to perform the 
prescribed operation as accurately as possible without 
resorting to the cut and try approach. 

1 . 1 Dielectric Resonator ; An Introduction 

A dielectric resonator is made of a low loss, high 
dielectric constant, and temperature stable ceramic material 
in a regular geometric forms like pillbox, disc, tubular, 
spherical etc. (Fig. 1.1) [12,14]. For a practical resonator 
the cylindrical (disc type) structure is preferred over 
other type of structures because cylindrical resonators do 
not have a large number of degenerate modes and are easier 
to fabricate and mount in microwave circuits [43] . In this 
thesis cylindrical dielectric resonator have been dealt in 
detail, due to its simplicity and the above mentioned 
advantages. Application of spherical DR, for miniaturization 
of DR filter for mobile communication has also been 
reported, recently [38] . 

In 1939 R.D. Richtmayer observed that unmatelized 
dielectric objects can function as electric resonators. The 
first exploratory activities on dielectric resonators 
occurred in the 1960s with the analysis of resonant 
frequency and modes, studies of dielectric resonator (DR) and 
circuit properties, and other related topics about coupling 
with external circuits. At that time, DR had been used in 
paramagnetic spin resonance experiments, but the lack of 
suitable material precluded its practical applications in 
microwave components. These early studies however helped in 
the understanding of DRs and stimulated the efforts to 
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develop suitable dielectric materials having desired 
properties for use in microwave applications such as 

a. high dielectric constant ed 

b. high unloaded quality factor, Qo 

c. low temperature coefficient of resonant frequency tf 

The value of Qo limited by major imperfections of the 
conductor and the dielectric losses occurring inside the DR 
material . The conductor losses are incurred due to small but 
non zero conductivity of the dielectric material . The 
damping of oscillations of electrons caused by the 

alternating polarization of dielectric material exposed to 
the time harmonic electric field at microwave frequency, 
results in the dielectric losses. Typical values of Qo for a 
DR is ~ 10,000. 

Temperature stability depends upon the effects of 

temperature on material and electrical properties of the 
resonator. These effects include the change in dimensions of 
the DR and the variation of dielectric constant ed. The 

resonant frequency of a DR mainly depends upon its 
dimensions and dielectric constant. These result in the 

deviation of the resonant frequency with temperature change. 
The sensitivity of resonant frequency with temperature is 
denoted by tf, which is also known as temperature 
coefficient of the resonant frequency. 

Various new compositions having acceptable values of the 
ed, tf, and Qo were explored with the advancement in 
material technology. These include ceramic mixtures 
containing Ti02, various titanates, zirconates and glass 
ceramic mixtures [12,13]. Temperature compensated composite 
structure resonators using ferrite and feroelectrics were 
also reported. This rapid development in material technology 
removed the restriction on the frequency range of the DR, 
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and now they can even be applied to microwave circuits 
operating in millimeter wave frequency range [13 , 35] . Table 
1.1 lists some important properties of commonly used 
materials for DR. 


Table 1.1; Material properties of some DRs [7,12] 


Material 

composition 

Qo 

ed 

Frequency 
range (GHz) 

Temperature co 
Tf (ppm/°C)“ 

{ZrSn)Ti 

Oxide 

>10,000 
(4.5 GHz) 

35.7- 

36.4 

1.45-8.86 

-3 to +9 

BaLnTi 

Oxide 

>3,000 
(3.0 GHz) 

80.0 

0.7-3.62 

-3 to +9 

BaZnTaTi 

Oxide 

>10,000 

(lO.OGHz) 

27.6- 

30.6 

5.5 - 
32.2 

-4 to +4 

Ba Titanium 
Oxide 

>6,000 

(4.5GHz) 

36.6- 

38.3 

0.8-5.21 

0 to 4 


1.2 Methods of analysis 

After availability of dielectric resonator with 
excellent properties, researchers started investigating the 
fields and the resonant frequency of DR in different circuit 
environments, using various approximate and rigorous 

analysis methods. 

1.2.1 Approximate Methods 

An accurate mathematical description of the 
electromagnetic field in a dielectric resonator is 
considerably more complicated than the field distribution in 
a hollow metallic waveguide resonator. Approximate methods 
treated DR as the dual of a metallic cavity. DR was modeled 
as a section of dielectric rod waveguide [14] . This model 
was simplified by introducing some approximation such as by 
assuming perfect magnetic conductor (PMC) wall over the 
curved surface of the DR. Cohn's method falls in the same 
class [14] . He modeled a cylindrical DR as a section of 
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cylindrical dielectric waveguide with air dielectric 
boundary in z (axial) direction (Fig. 1.2), then calculated 
resonance frequency and the coupling coefficient. The air 
filled sections on both sides of DR remain below cutoff due 
to low dielectric constant than that in the dielectric 
region. Thus the modes in the air filled regions are 
evanescent and fields decay exponentially in the axial 
direction away from each air-dielectric boundary. This model 
gives an error of more than 6% in resonant frequency 
calculation, but it helped in designing waveguide dielectric 
resonator filter. This model does not take into account of 
the field leaking through this surface, because it puts a 
fictitious open circuit (i.e. PMC) over the cylindrical 
surface of the DR , In an actual DR magnetic field is 
always present in the vicinity of the cylindrical surface 
containing the circumference of the DR due to air-dielectric 
boundary. Outside of the resonator magnetic field decays 
very rapidly with distance. A more accurate method, based on 
the variational procedure for computing the resonant 
frequency in TEoiS mode was later reported by Konishi et 
al. [16] . Iveland [15] computed resonant frequency for a 
square pillbox shape dielectric resonator using a similar 
procedure. A more realistic model for a DR placed under 
shielded condition was used by Itoh and Rudokos [17] . They 
made this model more realistic by removing the assumption of 
PMC of the Cohn's model and took into account the field 
distribution outside the resonator also. Even this method 
gives resonant frequency which is high by about 5% for 
isolated resonators. However, due to its simplicity in 
calculation it was adopted by many researchers. 

1.2.2 Rigorous Methods 

The approximate methods, described in the last section, 
are useful as aids in understanding the principles of 
operation of DRs, and they may be helpful to the initial 
stages of the analysis or design process, but, they fail to 
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provide the high accuracy usually required in present-day 
microwave circuit design. In an attempt to obtain accurate 
solution of the problem, DR placed in MIC environment was 
analyzed using rigorous analysis methods. Mode matching 
technique is used in some of these methods like Radial mode 
matching method and Axial mode matching method. Zaki and 
Chen [37] have applied the axial mode matching method to 
determine the TE, TM and hybrid modes of a DR placed 
symmetrically in a cylindrical cavity. Tsuji et al. [25] 
analyzed an isolated DR by expanding the fields inside and 

outside a pillbox resonator in terms of spherical waye_ 

fimctions and by matching the tangential field components of 
the two regions at the resonator's surface. Van Bladel 
applied perturbational approach by expanding the fields in 
asymptotic series in inverse power of the Tea. This method 
is exact only in the limit as ea goes to infinite. Glission 
et al.[24] used integral equation method by choosing a 
suitable Green's function which also constitute the kernel 
of the integral equation. The success of this method depends 
upon the choice of the Green's function which is simpler in 
the case of an isolated dielectric resonator where 
free-space Green's function is required. in case of shielded 
resonators, the determination of the Green's function is 
more difficult, particularly for more complex resonators. 
Glission simplified this task by subdividing the resonator 
into two regular partial regions and by erecting a electric 
or magnetic wall at boundary between the regions. For 
maintaining the correct field in both regions magnetic or 
electric current sheets were postulated on this wall . By 
enforcing the continuity conditions between the partial 
regions, an integral equation was obtained, which was solved 
using method of moments. The zeros of the resultant matrix 
determinant was calculated to obtain the resonant frequency. 

All these rigorous analytical methods are versatile and 
allow a more general range of applications. However, these 
methods are not used in practical design, because they are 
quite complicated as well computationally intensive. A 


y 
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method which offers both simplicity and accuracy is required 
for computer aided design purpose. A method using effective 
dielectric constant (EDO technique with DWM method is 
suggested by Mongia [20] for cylindrical DR. This method is 
easy to implement and gives results, which are accurate to 
1 % . 

1.3 Organization of the Thesis 

Design of microwave circuits employing dielectric resonator, 
like filters, oscillators, radar detectors, speed guns, 
cellular telephones etc., requires an accurate knowledge of 
resonant frequency of DR and its energy distribution in all 
the regions of the environment and its coupling with other 
circuit elements (transmission lines, other resonators 
etc.) . In this thesis these design aspects have been studied 
for a cylindrical dielectric resonator placed in 
Suspended- sxibstrate environment. Stress is over dominant 
TEoiS mode of the resonator. 

The thesis is divided into five chapters. Chapter 2 
discusses the use of EDC method combined with DWM method for 
calculating resonant frequency and field distribution of a 
cylindrical resonator for TEoiS mode. Chapter 3 presents the 
electrical and magnetic energies distribution in all the 
regions of the DR's environment. In chapter 4 coupling 
between two identical dielectric resonators is presented for 
three very important configurations, along with the 
evaluation of coupling coefficient between them. And the 
last chapter deals with the characterization of coupling 
between cylindrical resonator and two very important 
transmission line configurations like suspended stripline 
and inverted microstripline with finite thickness conductor. 
Following this chapter comes the conclusion and future scope 
these studies. 
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(a) Cylindrical 


(b) Rectangular box 


(c) Hollow cylinder 


^^9 1.1: Dielectric resonators in different geometrical shape 


Dielectric Resonator with Magnetic Wall 



Air filled Waveguide in Evanescent mode 


Fig 1.2: 


Dielectric resonator 
(Chon's model) [14] 


3s a section of Dielectric rod 




Chapter 2 


CALCULATION OF RESONANT FREQUENCY USING EFFECTIVE DIELECTRIC 

CONSTANT TECHNIQUE 

In the design of microwave circuits having DR, the first 
step is the analysis for finding the field distribution and 
resonant frequency of the dielectric resonator (DR) . This 
analysis forms the sxibject matter of the present chapter. 

Initially, the researchers tackled this problem by the 
approximate methods such as second order magnetic wall 
model [14] for finding simple design solutions for the lower 
order modes. Though approximate models are very efficient 
yet they are not accurate enough. This deficiency can be 
overcome by rigorous methods such as Method of moments and 
Green's function approach, but, at the price of increased 
complexity and computational time. Researchers are always in 
search of an efficient method capable of producing results 
comparable in accuracy to those obtained from rigorous 
methods, for most frequently used TEoi5 mode. 

To analyze the dielectric resonator in TEoiS mode, the 
use of Effective Dielectric Constant (EDC) technique applied 
with Dielectric Waveguide Model (DWM) method, is described 
in this chapter 

2.1 DWM method with EDC Approach 

2.1.1 DWM Method 

In DWM method a dielectric resonator is considered as the 
dual of a metallic cavity. Our Dielectric waveguide model 
method is similar to Itoh's method [17] except in field 
expressions selected for axial magnetic field (Hz) . The 
dielectric rod waveguide when truncated becomes a dielectric 
resonator, in same way as a truncated hollow waveguide 
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becomes a resonant cavity. For circular symmetric field 
configurations (lower frequency modes) , suitable field 
expressions can be found for axial field components (Ez and 
Hz) from the knowledge of field analysis for dielectric rod 
waveguide. When we apply the proper boundary conditions 
these field expressions give characteristic equation, the 
solution of which gives the resonant frequency. Analysis 
presented below is for TEoiS mode but it equally applies to 
TMoiS mode by virtue of the principle of duality. 

Maxwell ' s equations for TEon case 

For a source free, isotropic and homogeneous medium, the 
fundamental Maxwell's equations can be written as follows: 


VxE=- jwuH 


VxH= j o)c E 


( 2 . 1 ) 


and, 

V.E=0 

V.H=0 (2.2) 

Combination of above two equations gives us well known 
vector Helmholtz's equations given below ; 

V^E + k^E = 0 

V^H + k^H = 0 (2.3) 

where k is called the wave number of the medium and in 
general it is a complex quantity. For a dielectric medium of 
relative dielectric constant cr, the wave number k is given 
by: 

k^ = kf er (2.4) 
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where, ko = 2n fo^ noCo and fo is resonant frequency. 

To obtain the fields in all the regions of a structure, 
it is convenient to start from the axial component of the 
electric or magnetic field. Other field components can be 
obtained from this axial field with the help of Maxwell's 
equations . 

For TE^ mode (i.e. Ez = 0), We can start from the axial 
magnetic field i.e. Hz. Then vector Helmholtz's equation 
(2.3), reduces to scalar form as follows: 

V^Hz + k^Hz = 0 (2.5) 


we can solve this equation by well known method of 
separation of variables. The field components for 
TE^circular symmetrical case (d/dz = 0) are expressed as: 

T, TT a^Hz 

Er=0 Hr= 

V^dzdX 



E<p=Q 


Ez=0 Hz=cho sen. (2.6) 

Here, v is defined as the radial wave number, and in 
general is given by : 

2 2 2 

V = ±{k - kz ); where kz is the axial wave number. (2.7) 
Field distributions in the resonator 

The cross-sectional view of cylindrical resonator placed 
in suspended substrate environment is shown in Fig. 2.1. The 
resonator material is assumed to be a perfect dielectric 
characterized by the real scalar relative permittivity cd 
and the real scalar permeability ii (=fJ.o) . To simplify the 
analysis, this composite structure is divided into eight 
regions with a different set of field components. Here, the 
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proper choice of Hz in different regions depends on the 
requirements imposed by physical geometry and the 
analyticity of the field at the origin. It also should 
satisfy the wave equation (2.5). 

Inside the resonator, Hz should represent standing wave 
in z direction and wave number v = k.r should be real to 
fulfill the requirement of DR as a energy storing element. 
The solution of eqs.(2.5) in cylindrical coordinates, using 
method of separation of variables with above mentioned 
conditions gives the following form of expressions for Hz 
in region no . 1 

Hzi = {cos(|3iz) + Bi sin (|3iz) } Jo (krr) (2.8) 

where Jo is the Bessel function of the first kind and zero 
order and kr is related to pi by separation relation 

kr^ + pi^ = edko^ (2.9) 

Outside the resonator, Hz should be axially symmetric, 
but, decaying in nature in the direction away from the 
resonator for confirming the maximum energy storage ins\de 
the resonator. This decay in radial direction can be 
represented by modified bessel function, so, outside the 
resonator (in the radial direction) radial wave number 
becomes imaginary. Finally, normal magnetic field Hz should 
vanish at metallic surfaces. Again solving eqs(2.5) with 
these conditions, Hz in the rest of the regions (2 , 3 , 4 , 5 , 6, 7 
and 8) can be written as ; 

Hz2=B2.exp(-aadi) .sinh{aa (di-z) }Jo (krr) 

Hz 3=B3. {cos (|32Z) + B4 . sin 02Z) } Jo (krr) 

Hz4=B5.exp (-asdz) .sinh{as (dz+z) }Jo (krr) 

Hz 5=B6. exp ( -ttadi ) . sinh{aa (di-z) }Ko (kar) 
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Hz6 = B7.{cosOiz) + Bi . sinOlz) }Ko (kar) 


Hz7 = Be. {cos {|32z) + B4 .sinOzz) }Ko (kar) 

Hzs = Bg.exp (-asdz) . sinh{as (d 2 +z) }Ko (kar) (2.10) 

where Ko is the modified Bessel function of the second 
kind and zero order, and separation relations in these 
regions are 


1 2 
ka 

= 

^ 2 
- ko 

(2. 

.11) 

1 2 
ka 

= p2^ 

1 2 

“ €sl . ko 

(2. 

.12) 

1 2 
kr 

- aa^ 

II 

o 

to 

(2. 

.13) 

1 2 
kr 

- OL^ 

1 2 

= Cs2.ko 

(2. 

.14) 


The transverse field components in various regions are 
found by substituting (2.8) and (2.10) into (2.6), and are 
given as 
Region 1 

E 01 = _ (cosOiz) + B1 . sin Oiz) } Ji (krr) 

kr 

Hri = (j3i/kr) {sin(/3iz) - Bi .cos (|Siz) } Ji (krr) (2.15) 

Region 2 

E<p 2 =- - B 2 .exp (-aadi ) . sinh{aa (di -z) } Ji (krr) 

kr 

Hr 2 =B 2 . (aa/kr) .exp (-ttadl ) .COShftta (dl -z) } J 1 (krr) 

(2.16) 
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Region 3 


E^3=- B3. {cos {/S 2 z) + B4 . sin (/32Z) } Ji (krr) 

kr 

Hr3=B3. 02 /kr) , (sinOzz) - B4 .COSO 2 Z) }jl (krr) (2.17) 

Region 4 

E^ 4 = _ B5.exp(-asd2) .sinh{as (d2+z) }Ji (krr) 

kr 

Hr4=-B5. (as/kr) .exp (-ttsdz) .cosh{as (d2+z) } Ji (krr) (2 . 18) 
Region 5 

E<^5 = B6.exp(-aadi) .sinh{aa (di-z) }Ki (kar) 

ka 

Hrs = -Be . (aa/ka) .exp(-aadi) .cosh{aa (di-z) }Ki (kar) (2.19) 

Region 6 

E06 = B7.{cosOiz) + Bi . sin Oiz) }Ki (kar) 

ka 

Hr6=-B7 . (jSi/ka) . {sin(piz) -Bi . COS (piz) }Ki (kar) (2.20) 

Region 7 

E07= B8.{cos02z) + B4 . sin (j32z) }Ki (kar) 

ka 

Hr7=-B8 . (|S 2 /ka) . {sin(/S2Z) -B4 .COS (/S2z) }Ki (kar) (2 .21) 

Region 8 

E^8=^ — — B9 .exp(-o:sd2) .sinh{as (d2+z) }Ki (kar) 

ka 

Hr8=B9 . (as/ka) .exp (-ttsdz) .cosh{as (d2+z) }Ki (kar) (2.22) 

The value of unknown coefficients Bis, can be evaluated 
by applying the boundary conditions at the interfaces of two 
regions, and are given in appendix A. 
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Resonant frequency £o 


If kr and pi are known the resonant frequency fo can be 
calculated from eqs, (2.9). The concept of dielectric 
waveguide model (DWM) is used for determination of kr and pi 
[20] . The kr and pi are determined as follows. 

Determination of kr 

The resonator is first assumed to extend to infinity in z 
direction (Fig. 2.2(a)) for determining kr. We assume that 
this infinite dielectric cylinder carries TEoi as the 
propagating mode, and that kr is the radial wave number for 
this cylinder for r < d. Applying the continuity condition 
for tangential electric field and normal magnetic field 
components at the boundaries of region 1, gives following 
characteristic equation: 

.. + K^lkad) ^ Q (2 23) 

krjo (krd) kaKo(kad) 


Radial wave number kr is related to ka by 

ka = { (Cd - 1.0) ko^ - kr^}^^^ (2.24) 

kr and ka can be computed by solving Eqs. (2.23) and (2.24) . 
Determination of pi 

The structure of Fig. 2.2 is assumed to be a radial triple 
slab guide (Fig. 2.2(b)), along which the TEoi guide mode is 
propagating. In other words it is assumed that the axial 
wave number pi is the same as that of a triple layer slab 
guide structure, obtained by extending the resonator to 
infinity in the radial direction. When we apply the boundary 
conditions at the interfaces of these regions we get one 
more transcendental equation, given by 
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^ 2 [M 2 {tan (|32h/2) - tan02d3)} + tan 02 h/ 2 ) .tan02d3) +1.0 ] 


[M 2 + tan(/32d3) - tan(|S2h/2) + M 2 tan 02h/2) tan{/32d3) ] 

/Si [tan^ (/Sih/2) + 2 Mi tan(/3ih/2) - 1.0] 

= 0 (2.25) 

[Ml - Mitan^ (^ih/2) + 2 . OtanOih/2) ] 


where , 

Ml = (^i/aa) . tanh (aaha) ;and 

M2 = (/S2/as) .tanh(ashs2) (2.26) 

We get the value of kr by solving transcendental 
equation (2.23) in conjunction with equation (2.24) . Then we 
can solve the equation (2.25) for the height of the DR for 
specified resonant frequency fo, after substituting the 
values of |3i , /32, aa, and ots from separation relations 
(2.12) - (2.14) . 

For given dimensions of the DR, the resonant frequency 
of the DR is obtained by satisfying two transcendental 
equations (2.23) and (2.25) simultaneously, for suspended 
substrate as well as microstrip (csi = es2 = Cs) 
environment . 

2.1.2 Use of EDC Approach 

The DWM method gives higher value of the resonant frequency 
fo. In higher frequency range specially for millimeter waves 
this error is very large. Recently Mongia [20] suggested 
Effective Dielectric Constant (EDC) approach for calculating 
radial wave number kr, which greatly enhances the accuracy 
of DWM method. According to this approach radial wave number 
kr is assumed to be the wave number of an infinitely long 
cylindrical waveguide having the same radius as that of the 
DR, but having a dielectric constant Ceff which is less than 
that of the resonator. In this method kr and ka are related 
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by 


ka = {(eeff - 1.0) kl - k? (2.27) 

An approximate expressions for eeff [20] is given as follows 
Ceff = (eeff* + eeff*”) /2 (2.28) 

where 


a 

eeff 


cd - (ed-eeff ) . (hsi+hs 2 ) /d if (hsi+hs 2) < d 
eeff otherwise 

(2.29) 


= ed - (ed - eeff ) ?ha/d if ha < d 


eeff 


= Ceff 


otherwise 


(2.30) 


Ceff* is obtained by 


eeff = 


kr (ed) 

1 2 
ko 


(2.31) 


Here, kr(ed) is the same radial wave number given in eqs 
(2.24) . 

Thus the values of kr, and ka are modified by EDC method 
and also, the values of other eigen values are 
correspondingly modified. When (hsi+hsa) = ha = 0, the value 
of eeff reduces to ed and this method gives the exact value 
of resonant frequency for a dielectric post resonator. When 
(hsi+hs2) and ha both tends to infinity with csi = es2 = 
1 . 0 , eeff reduces to kr(cd)/ko, and gives the expression for 
an isolated resonator. 
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Results obtained by EDC approach are comparable to those 
obtained by rigorous methods . 

2.2 Numerical Implementation 

Above analysis is very general in nature, that is applicable 
for suspended sxibstrate environment as well as for MIC 

environment (By putting esi = cs 2 = cs and hsi = hs 2 = 
hs/2) . Based on this analysis software programs for both, 
pure DWM and DWM with EDC approach, named f-DWM.f and 

f-EDC.f, are developed for the calculation of the resonant 
frequency for a given dimensions of the DR, and the 

supporting structure. 

In these programs computation is started with finding 
the roots of equation (2.23) using Newton-Rapshon method 
with following initial approximate value for kr (ea) : 

kr(ed) = [0,951 Xoi+0.222{ (Cd-l) (kod) 0.951 Xoi}^^^]/d 

(2.32) 

where Xoi = 2.405, and is the first root of the equation 

Jo (x) = 0 . However calculation of kr in above closed form 

expression requires the value of fo also. A wise guess for 
this value reduces the number of iterations required for 

convergence of the solution of the equation (2,23) and 

consequently the computational time of the programs. This 
value is obtained from following empirical relation, 
obtained by curve fitting process 

fo(GHz) = 36.0 (^ + - + 1-^2^ — )/4~cd (2.33) 

here h is height in mm, d is radius in mm, and ea is 
dielectric constant of the DR, And ha is the distance in ram 
of the top conductor wall from the top of the DR. This 
relation gives the starting value of f o accurate to ± 4% . 

Other eigen values and the height of the resonator can 
be computed, when kr and ka are determined. Repeated 

iterative procedure is used to get the value of fo within 
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tolerable accuracy. 


2.3 Numerical Results 

2.3.1 Comparison with Experimental Data 

The results obtained for fo of a DR from above programs are 
presented in Table 2.1 along with the experimental results 
[20] in MIC environment, because of non availability of 
experimental data in the literature for suspended- sxibstrate 
environment. But, as mentioned earlier, it is a general 
analysis which is applicable for both types of the 
environment, so the verification in either case with 
experimental results validate the analysis. 

A comparison of theoretical and experimental results 
shows a good agreement in all cases, the difference between 
the two results is less than 1% in most cases. 


Table 2.1: Comparison of EDC method with Experimental 
Data [201 (d = 4.55 mm, ea = 37.1, ha infinite) 


hs 1 

Gs 1 


Resonant frequency 

fo (GHz) 


=hs 2 

=:CS 2 













=hs/2 

(mm) 

=C r s 

h = 2 

1 

. 0mm 

h = 4 

. 04mm 

h = 6. 

96tnm 




Exp. 

Theory 

Exp . 

Theory 


0.19 

2.22 

8.81 

8.89 

6.46 

6.46 

5.43 

5.37 

0.4 

2.22 

8.3 

8.26 

6.3 

6.24 

5.34 

5.28 

1.19 

2.22 

7.39 

7.44 

5.91 

5.87 

5.17 

5.09 

1.58 

2.22 

7.3 

7.3 

5.82 

5.79 

5.14 

5.05 

0.32 

9.6 

8.4 

8.43 

6.31 

6.31 

5.37 

5.31 
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Table2.2 Comparison of EDC and DWM methods with Rigorous 

methods 


Structural parameters Resonant frequency 

(GHz) 


ed 

es 1 = 

es2 

d 

(mm) 

h h 

(mm) 

s 1 =hLs 2 

(mm) 

ha 

(mm) 

Rigorous 

methods 

EDC pure DWM 

tech method 

38.0 

1.0 

2.86 

2.38 

0.36 

inf. 

9.83 

9.82 

10.22 

38.0 

1.0 

2.86 

2.38 

1.43 

inf . 

9.16 

9.14 

9.85 

35.2 

9.6 

2.43 

1.81 

0.32 

inf . 

12.4 

12.37 

12.87 

34.19 

9.6 

7.49 

7.48 

0.35 

0.72 

4.35 

4.33 

4.35 

34.21 

9.6 

7.00 

6.95 

0.35 

1.25 

4.52 

4.50 

4.53 - 

36.2 

1.0 

2.03 

5.15 

1.465 

2.93 

10.50 

10.36 

10.85 

36.2 

1.0 

4.00 

2.14 

2.215 

4.43 

7.76 

7.68 

8.365 


2.3.2 Comparison with DWM and Rigorous Methods 

Using this EDC method the results were computed for 
microstripline configuration, and the results are compared 
with approximate [17] and rigorous methods available in the 
standard literature. Table 2.2 shows this comparison. It is 
to be noted that the present technique gives results having 
an accuracy comparable to that of rigorous methods . The 
difference never exceeds more than 1%. In Fig. 2.3, results 
obtained using pure DWM method are also compared. It is seen 
that when ha is small as compared to d then the error of the 
DWM method in calculating fo is very small. This is expected 
since ceff for this case almost reduces to cd. 

2.3.3 Effect of Design parameters on f o 

Various electrical and physical parameters of the circuit, 
which are under the control of the designer at time of the 
fabrication the circuit, may change the resonant frequency 
of DR in the circuit. . Information presented in Figs. 2.4 
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to 2.6 shows the impact of these control variables over 
resonant frequency fo. 

Fig. 2.4 shows the effect of separation (ha) between top 
conductor plate and the top surface of the DR, over fo. 
Variation of resonant frequency with changing ha can be 
explained by the cavity perturbation theory. When a metal 
wall of a resonant cavity is moved inwards, the resonant 
frequency of the cavity will decrease if the stored energy 
of the displaced field is predominantly electric. Otherwise, 
when the stored energy close to metal wall is magnetic, as 
in the case of shielded TEoiS mode, the fo will increase 
[7]. Fig. 2.5 is showing the effect of the dielectric 
constant (ea) of DR, on its resonant frequency. We find that 
the resonant frequency decreases with increase in the value 
of Gd. Hence the material of DR should be temperature stable 
to get frequency stability. Fig. 2.6 shows the variation of 
fo with substrate thickness (hsi) for different values of 
substrate dielectric constant (esi) . 

Resonator Tuning : 

Tuning a dielectric resonator means adjusting its 
resonant frequency by mechanical or electrical means. 
Several techniques can be used to achieve this like 

# changing the thickness of the dielectric resonator. 

# perturbing the fringing fields outside the DR via 
screws, tuning stubs, or ferromagnetic material, etc. 

# perturbing the DR's fo by the so-called double 
resonator configuration, where two halves of the DR act as 
one. A tuning range (up to 20%) can be achieved with little 
degradation in Q factor by this method. 

# Optical tuning, magnetic tuning and electric tuning is 
also used. 

The choice of a tuning method requires the consideration 
of the following factors.* 

a. the mode to be tuned, 

b. the mechanical configuration of the final device, and 

LIBH Aft'S 

i s. T., KANPlfw 

~ r-Trttn 



c. the tuning range desired and its direction (lower or 
higher frequencies) 

For example, the tuning of TE modes is easily achieved 
by metallic or dielectric tuning stubs placed perpendicular 
to the dielectric resonator's top surface. Fig. 2.7 shows a 
schematic of the general effect of these methods. Metal 
tuning stubs (a) moving towards the DR will pull the frequency 
up and dielectric tuning stxibs (b) will push the frequency 
down. 

Approximately 20% of fo tuning can be achieved by these 
methods. It is, however, good practice to restrict this 
amount to below 5% to minimize degradation of both 
temperature coefficient and unloaded Q. 
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fo{GH2) 


d = 0.96482inm h = 0.68395inm 
e«i = 2.22 hsi = 0.79375mm 

~ 1.00 lis 2 = 2.00mm 

ed = 29.4 


Fig. 2.4: Effect of changing height of top conductor wall (ha) 
on resonant frequency fo 








Fig. 2.6: Effect of changing thickness and dielectric constant 
of the suspended substrate. 



Jlia — 2.7f Tiini no bv stub method [12] 






Chapter 3 


ENERGY DISTRIBUTION OF DIELECTRIC RESONATOR IN 
SUSPENDED SUBSTRATE ENVIRONMENT 

The relative energy distributions in different regions of 
the circuit structure involving dielectric resonator, is 
very useful in deciding how to couple to the resonator. This 
chapter will deal with electrical and magnetic energy 
distribution for TEoiS mode at resonance. 

3.1 Electrical Energy Distribution 

Total electrical energy stored in a structure given in 
Fig. 2.1 is the sum of individual electrical energy stored 
in each region. It may be calculated from field expressions 
derived in chapter 2 . 


8 

We = E Wei 
1 = 1 


(3.1) 


. th 


here Wei is the electric energy stored in the i“‘ region of 
the structure occupying volume V and having dielectric 
constant eri . Wei is given by 


Wei 


CoGri 


S |El| dv 
V 


(3.2) 


For TEoiS mode only E^ component of electric field exists so 
for different regions we get stored electrical energy as 
follows : 

For region 1 


Wei = Xj X^’^ lE^ij^r dr d(p dz 

Z 0 0 “ 11/2 


(3.3) 
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region 2 


We 2 = Xq X ^2 \EcP 2 \\ dr dz (3.4) 

region 3 

We3 = ~^Sl X:^^ |E^3|^r dr d<f> dz (3.5) 

region 4 

We4 = ££|l2j-d J.2n j.-d3 |e^4|2j. (Jj- d^ dZ (3.6) 

region 5 

Wes = — 1^ d^ dz (3.7) 

region 6 

We6 = \E<p6\^r dr <l<p dz (3.8) 

Z a o “n/2 ‘ * 


region 7 


We7 = Xl^^ |E<^7|^r dr d<t> dz (3.9) 

region 8 

Wes = J'Id 2 |E^8|^r dr d^ dz (3.10) 

here ts is the saturation distance, i.e. after this distance 
in radial direction energy is zero. 

After putting the values of fields from chapter 2, and 
solving "integrals we get the expression for these electrical 
energies as follows: 

Wei = AO.A1234,Cd.A16 (3.11) 
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(3.12) 


We2 = A0.A1234.B2. t same ( ha , dl , <Xa ) 

We3 = A0.A1234.Csi.B3.A37 (3.13) 


We4 = 

2 

: AO .A1234. CsZ.BS. tsame{hs2,d2,as) 

(3.14) 

WeS 

2 

= AO s A3 678 s 1^6 {ildi g dl f CX&) 

(3.15) 

We6 

= A0.A5678.B7.A16 

(3.16) 

We7 

= A0.A5678.Csi.B8.A37 

(3.17) 

We8 

= AO .A5678. Cs2.bI. tsame (hs2, d2, tts) 

(3.18) 

where , 



AO 

= n.Co. (CJo^to) ^ 

(3.19) 

A1234 = lej/kf 

(3.20) 

A5678 = lek/kl 

(3.21) 

Here constants Bis are the same as in 

chapter 2 , and are 

given in 

appendix A. Other constants are 

given in appendix B 


and C. 

3.1.1 Electric Energy filling factor 

Energy filling factor is an important quantity used in 
concern with the analysis of partially filled cavities. This 
is defined as the ratio of the electric energy for a given 
region to the total electric energy stored in the system. 
Based on the above analysis for electric energy, a program 
is developed to find the energy in all the regions. Results 
obtained are tabulated in terms of energy filling factor 
for each region for different dielectric constants of the 
resonator and substrate. 
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Table 3.1: Electric Energy Filling Factor for various 
regions of resonant structure in suspended substrate 
Environment 


(d=0 . 9648mm, 

£82=1 .0) 

h=0 . 684mm, 

hsi=0 

. 79mm, 

hs2=2 . Omm, 

lla=3 . 

Omm, 

Cd 

Cst 

Energy filling factor (= 

12 3 4 

Wei/We)% for 

5 6 

region 

7 

no. 

8 

29.4 

2.22 

97.5 0.44 

1.04 

.005 

.096 .733 

0.23 

.001 

35.0 

2.22 

97.9 0.36 

0.86 

.004 

0.08 0.6l 

0.19 

.008 

29.4 

10.0 

96.1 0.42 

2.49 

.004 

0.08 .495 

0.45 

.000 


It can be observed from Table 3.1. that for ea = 35.0, most 
of the electrical energy (97.8%) is confined within the 
resonator. Energy leakage in region increases when 
dielectric constant of the DR decreases. Similarly using 
higher dielectric constant substrate, energy contents in the 
region 3 and 7 increase. But, in all cases more than 96% 
electrical energy is stored in the resonator. 

3.2 Magnetic energy distribution 

As we have seen from previous description that for high 
dielectric constant resonators, more than 96% electrical 
energy is confined within the resonator for TEoiS mode, 
hence, electrical coupling of the DR with other circuit 
components is very week. We can also find magnetic energy of 
the structure in the same way as we have done for electrical 
energy. The net magnetic energy is the sum of magnetic 
energies of individual regions. 

8 

Wm = E Wml (3.34) 

1 = 1 

Here Wmi is the magnetic energy stored in the i*^*' region of 
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the structure {Fig. 2.1), occupying volume V, and having 
relative permeability iin . Wmi is given by 

Wmi = s |Hi|^ dv (3.35) 

V 

For TEoiS mode only Hr and Hz components of magnetic field 
exists. So, 

|Hi|^ = Hi . Hi = nfi + nil (3.36) 

And for different regions we can write as follows: 

For region 1 

Wmi = {H?i+ nil} r dr d^ dz 

2 0 0 -h/2'- ■' ^ 

region 2 

Wm 2 = {H?2+ rIz} r dr d<p dz 

2 0 0 h/2 

region 3 

Wm3 = Hi3} r dr d^ dz 

2 0 O -d3 ^ 

region 4 

Wm4 = X"**^{I^4+ 1 ^ 4 } r dr d<^ dz 

2 0 0 -d2 ** ^ 

region 5 

Wms = -^X^'X^^ Xf„ {H?5+ hIs} r dr d^ dz 

Z d 0 h/2 

region 6 

Wm6 = -i^X^®X^^ {R^6+ hIs} r dr d(p dz 

2 d u —h/2 
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region 7 


Wm7 = — ^ ^2 

region 8 

Wm8 = ^2 ^ 

Here, ts is the saturation distance, i.e. after this 
distance in radial direction energy is zero. 

After putting the values of fields from chapter 2, and 
solving integrals we get the expression for these electrical 
energies as follows: 

Wmi = A0{A1234.A16+ pi .A1234. A16} 

^ 2 2 
Wm2=A0. B2{A1234.t same (ha, dl , OCa) +OCa . A1234 . t same (ha, dl , Cta) } 


Wm3 = A0.B3{A1234.A37 + p2.A1234.A37} 

Wm4=^0. 35 {A1234 , tsame (hs2, d2, OCs) +OCb .A1234* t same (hs2, d2, CCs) } 

Wro5=J[0, B6 {A5678 . tsame (ha, dl , OCa) +C)Ca . A5678 * t same (ha, dl | CCa) } 

Wm6 = A0.B7{A5678 .v115 + .A567S . A15} 

Wm7 = AO. bI {>15678. A37 + /sl . A5678 . A37} 

A 2 2 * *1 

Win8=A0. B9{A5678 . tsame (hs2, d2, OCs) +tts .A5678 . t same (hs2, dz, Ots) } 

(3.37) 

where 

AO = jLioH 
AI234 = Ijo/k> 

A5678 = Iko/ki (3.38) 

and other constants are given in appendix B and C. 
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3.2.1 Magnetic Energy Filling Factor 


This factor for a region represents the part of total 
magnetic energy in that region. For TEoiS mode it is of 
great importance, because, the coupling is mainly magnetic 
energy coupling in this mode. Table 3.2 is showing the 
magnetic energy distribution in all the regions of the 
stmacture . 


Table 3.2: Magnetic Energy Filling Factor for various 
regions of resonant structure in suspended substrate 
Environment 


{d=0 . 9648mm, 

es2=l .0) 

, h=0. 684mm, 

hsi=0 

. 79mm, 

hs2= 

=2 . 0mm, 

ha=3. 

0mm, 

ed es 1 

Energy filling factor (= 

12 3 4 

Wmi/Wm)% for 

5 6 

region 

7 

no. 

8 

29.4 2.22 

58.1 

13.78 

14.2 

.15 


8.67 

2.55 

.03 

35.0 2.22 

58.3 

13.76 

14.1 

.14 

2.47 

8.7 

2.5 

.03 

29.4 10.0 

63.8 

15.81 

7.04 

.17 

2.55 

9.54 

1.1 

.03 

We can 

observe 

from 

above 

Table 

3.2 

that in any 

case 


near about 38% of the magnetic energy is distributed out of 
the resonator. Therefore, strong coupling between resonators 
and other components can be achieved by placing the 
component in the region of higher magnetic energy. For this 
purpose region no. 2 and 6 can be used for strong coupling 
and for week coupling requirements region no. 4 or 8 can be 
used. 

The information presented in this chapter is very useful 
from design point of view. It will be used in proceeding 
chapters . 
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Chapter 4 


COUPLING BETWEEN TWO SYMMETRICAL DIELECTRIC RESONATORS 

IN SUSPENDED SUBSTRATE 


The design of microwave circuits, involving dielectric 
resonators like band pass or band reject filters, basically 
involves two design issues: the first one is the 
input/output coupling and the second is the interstate 
coupling. Interstate coupling, which provides transfer of 
energy from one resonator to . another resonator either 
directly or through another transmission line, is presented 
in this chapter. 

Analysis presented in the next chapter may be extended 
for calculating the coupling between two resonators coupled 
via a section of suspended stripline, micro stripline or 
inverted microstrip line. In this chapter we derive the 
analytical formulation for the coupling coefficient between 
two directly coupled DRs, for three very important 
configurations viz. edge coupled, off-layered and broad-side 
coupled resonators. 

4.1 Direct coupled Cylindrical Resonators 

Two cylindrical resonators may be coupled directly in 
TEoiS mode by placing them over a suspended substrate 
adj acent to each other as shown in the Fig . 4.1. The degree 
of coupling is mainly determined by the lateral distance S 
between them. As both resonators are assumed to be operate 
in the TEoi5 mode, interaction between them occurs through 
magnetic field lines, because, most of the electric field is 
confined within the resonators itself. Due to this magnetic 
interaction, resonators may be assimilated as magnetic 
dipole, formed by small current loops. 
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Fig, 4. 2 shows the side view of two cylindrical resonators 
placed at a distance of S. From this figure we observe that 
first resonator is placed in the sixth region (refer to Fig. 
2,1) of the second resonator and vice versa. As we have 
Assumed both resonators to be identical (in dimension and 
material properties) , each of them may be represented by a 
conducting loop having self inductance Lr resonating with 
capacitance Cr at the resonant frequency fo. Magnetic field 
interaction between the resonators is characterized by 
mutual inductance Lm. In circuit theory, the coupling 
coefficient k is defined as 


k 


Lm 

Lr 


(4.1) 


Here, we are ^ representing the resonators by equivalent 
current loops and the voltage induced in loop 2 is due to 
current In flowing in loop 1 . In terms of circuit 
quantities this induced voltage may be expressed as 

V2 = - Lm = -jwoI«.In (4.2) 

Value of the magnetic energy stored of the resonator 1 is 
given by 


Wml = Lrin (4.3) 

From eqs. (4.1) and (4.2) we get 

Lm _ _ V2lrl f. 

Lr “ 2jWoWml 

Induced voltage in loop 2, according to Faraday's law, may 
be written in terms of field quantities as 

V2 = - i E^dl = juoiioS H^ds (4.5) 

C 

here, E^and H^are fields in place of the second resonator 
due to first resonator, and surface 's' is the area in 
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second DR and substrate through which magnetic lines crosses 
perpendicularly. ' c.' is the contour enclosing the surface 
's'. Here it is assumed that voltage induced in the one 
resonator is only due to the current flowing in the other 
resonator. 

From reaction concept it comes that 

Inf^E^dl = dv (4.6) 

where jl is the current density in the volume of second 
resonator. 

From eqs (4.1) to (4.6), we get the expression for coupling 
coefficient as follows: 

XXX dv 

k = — (4.7) 

2j 0)0 Wml 

If current loops are replaced by original DRs in the above 
analysis then, current density J becomes equal to 
polarization current ^ present in volume v. In the present 
context jIp is expressed as 

Jp = jwoeo{ (ed-l)E^ +(esi-l)E^ +{cs2-l)E^} (4.8) 

where E^ is electric field of the second resonator in volume 
V. For TEoiS mode only E<p component of electric field exist, 
so, eqs (4.8) reduces to 

jIp = jt<)oeo{ (cd-l) e| +(esi-l)E^ +(es2-l)E| }l^ (4.9) 

For calculation of coupling coefficient using (4.7) at given 
resonant frequency fo, we need to find out two things. They 
are the magnetic energy stored (Wmi) due to excited 
resonator and the volume integral which changes according to 
placement of second resonator. Wmi which reman constant for 
fixed structural parameters, can be obtained from previous 
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chapter. Here we will present the volume integral for three 
stated configurations. 

4.1.1 Ed ge -coupled Resonators 

For edge coupled resonators (Fig. 4.2), the volume integral 
can be written in the form of 

SSSE^^ dv = [icd-l)B7A16 + (Csi-1)B3B8J137 

+ (es2-l)B5B9 tsame (hs2,d2,as) ] Ijk 

(4.10) 

The value of constants are given in Appendix A, B and C. 

By sxibstituting the values of volume integral from 
egs.(4.l0) to eqs, (4.7), we get formulation for coupling 
coefficient between two DRs placed on the upside of the 
suspended substrate (Fig, 4.1). 

4.1.2 Coupl ing between two DRs placed on Opposite sides 
of Suspended substrate 

We have extended this theory for getting coupling 
coefficient between two DRs, placed either side of the 
suspended substrate (Fig. 4.4) by replacing the two DRs with 
two equivalent current loops. The current flowing in one 
loop will generate the magnetic flux, which will be linked 
to the second loop through the substrate. The change in 
magnetic flux will give rise to induced voltage in second 
loop. We have assumed both DRs to be identical in properties 
and dimensions, and the formulation is derived for 
symmetrical location of DRs on the suspended substrate 
(i.e., hs 2 = ha + hsi, and es 2 = 1.0), for two 

configurations of Fig. 4.4. The volume integral for above 
mentioned configurations can be given as 

XJ'Xe\j dv = (cd - 1.0) .e'“-^^.B9.Als 

+ (ers - 1.0) .B3.B8.A37d ] .Ijk (4.11) 

for off-layered case. 
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and, 


SSS^-H dv=2Tr 


jCJoCo (tJoHo) ^ 
kr .kr 

+ (es 


[ (ed-i.O) 

1 - 1.0) .B3. A37(f] .le J 


(4.12) 


for broad side case. 

The value of constants appearing in these expressions are 
given in appendix A, B and C. 

By substituting the value of volume integrals from eqs 
(4.12) and (4.13) into the eqs (4.7), we can obtain the 
coupling coefficient for two cases given in Fig. 4. 4. 


4.2 Numerical Results 


We have validated the present method, by calculating the 
coupling coefficient between DRs for given structural 
parameters [30] , in microstrip line as a special case of the 
suspended stripline. In our present calculation we 
calculated frequency and other eigen values from EDC method. 
These theoretical results shows a good agreement with the 
experimental data presented in [30] and can be seen in Fig. 
4.3. 

Fig. 4.5 shows the variation of the coupling coefficient 
between two edge coupled dielectric resonators placed in 
suspended substrate, as a function of edge to edge 
separation 'S' between them, for Ka band. It is clear from 
this graph that the coupling coefficient k is very strong 
function of separation S between the resonators. Hence, 
spacing between resonators should be determined very 
accurately in practical applications. 

Fig. 4.6 is showing the coupling coefficient variation 
with change in separation (S) between two Drs, for off 
layered coupling (Fig. 4.4a) . We see that the coupling 
coefficient k is a strong function of the separation between 
the DRs, and it decreases very rapidly as we increase the 
separation between two DRs. It also shows that for same 
dimensions of the structure, the coupling increases as we 
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increase the dielectric constant (ers) of the substrate. We 
also find that the coupling is week even for zero 
separation ( S=0 ) , compared to that in edge coupled case 
(Fig. 4. 2) . 

Fig. 4.7 is showing the change in the coupling coefficient 
with change in the width of the substrate (hsi) , for the 
Broadside coupling (Fig. 4. 4b). It shows that the coupling is 
decreasing with increase in the thickness of the substrate. 
It also shows that very tight coupling can be achieved in 
the Broadside case. 

We conclude from the above figure that off layered 
coupling between DRs (Fig. 4.4a), can be used for the 
purpose of week coupling and the Broadside stmicture (Fig. 
4.4b) can be used for very tight coupling. The most 
important advantage of these structures is that we can 
miniaturize the microwave circuits, involving the DRs for 
tight as well as week coupling purposes. 
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Fig. 4.1; Two directly coupled resonators on the same side of 
suspended substrate 



Fig. 4.2: Side view of two edge coupled dielectric resonators 






Fig. 4.3: Comparison of theoretical results with experimental 
data for microstrip configuration 
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Chapter 5 


COUPLING BETWEEN DIELECTRIC RESONATOR AND STRIP-CONDUCTOR 
WITH FINITE THICKNESS ON SUSPENDED SUBSTRATE 


As we have mentioned in earlier chapter that there are two types of 
viz. interstate' 

f 

coupling* and input /output coupling . The 

input/output coupling is the coupling between the resonator 
and a transmission line which either excites the resonator 
or extracts the electromagnetic energy from the resonator. 

Many types of transmission lines like co-axial probe, 
current loop, waveguides, stripline, microstrip line, etc. 
can be used for this purpose. Coupling between a cylindrical 
dielectric resonator (DR) and a strip conductor with finite 
thickness (t) , for TEoiS mode in suspended substrate 
environment, is the subject matter in this chapter. 

5.1 Methods of coupling for TEoi5 and TMoiS modes 

The TEoiS and TM 016 are most commonly used modes, which can 
be easily excited in dielectric resonator using microstrip 
line, fin line, magnetic loop, metallic and dielectric 
waveguides, microstrip line, suspended stripline etc. 

[ 8 ]. Field configurations in these modes are such that 
resonator operating in any of these modes can be modeled in 
terms of low frequency lumped elements. 

Fig. 5.1(a) shows the magnetic field distribution for 
TEoiS mode and various methods for coupling of this mode. 
Electric field lines for this mode are simple circles 
concentric with the axis of the cylindrical resonator . When 
the dielectric constant is around 40, more than 96% of the 
stored electric energy as well as a great part of magnetic 
energy (typically 60%) is confined within the resonator. The : 

remaining energy, distributed in the surrounding of the DR, 
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is decaying rapidly with distance away from the resonator 
surface. As seen from this figure, for a distant observer 
this mode appears as a small magnetic dipole formed by a 
small coil, and for this reason some authors call it a 
magnetic dipole mode instead of using the term TEoiS 
introduced by Cohn [14] . 

Coupling to this mode is often accomplished through the 
magnetic field via a small horizontal loop placed in the 
plane parallel to the flat surface of the DR or by placing 
the end face of the DR near a strip conductor of planner 
transmission line like microstrip line or suspended 
stripline (Fig. 5.19(a)), so that the magnetic lines of the 
DR link with those of the loop or the strip-conductor. 
Coupling of a resonator to a waveguide operating in the TEoi 
mode can be accomplished by placing the resonator axially 
within the waveguide. Coupling to this mode via electric 
field can also be obtained using a small horizontal dipole 
or a bent mono pole as illustrated in the figure. 

TM 016 mode is dual of TEoiS mode. Role of electric and 
magnetic field gets interchanged in this mode (Fig. 5.1(b)). 
Magnetic field is well confined within the resonator and 
electric field is strong near the axis of the resonator. 
Coupling to this mode can be accomplished in different ways 
as shown in Fig. 5. 1(b). A probe, well aligned with the axis 
of the resonator effectively couples the electric field. 
Coupling to strip conductor of planner transmission line can 
be accomplished by placing the resonator on its side with 
its axis of rotation parallel to the strip conductor so that 
the magnetic field lines loop around the strip conductor. 
The resonator may be coupled to the TEoi waveguide mode by 
placing it on its end in the waveguide as implied by the 
orientation of the waveguide in the figure. 


5.2 Coupling between DR and planner Transmission lines 

Mlcrostrlp line is the most popular transmission line used 
with cylindrical resonators, in microwave integrated 
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circuits (MICs) , operating at lower frequency band. But, at 
higher frequencies, specially in millimeter (mm) wave range, 
the microstrip line fails to satisfy the requirements of 
losses, fabrication tolerances and handling fragility etc., 
and also, it becomes very dispersive in nature. The 
supplements to these requirements are suspended stripline 
and microstrip line, for mm wave applications. Hence it is 
necessary to characterize quantitatively this coupling 
between a DR and two important transmission line structures 
like suspended stripline and inverted microstrip line. Our 
next sections are devoted to this transmission line to 
dielectric resonator transition. 

One way to incorporate the dielectric resonator in a 
microwave circuit is by placing it on the suspended 
substrate having the strip conductor of a planner 
transmission line, as shown in Fig. 5.2. 

When a DR is placed in the vicinity of a transmission 
line, as shown in Fig. 5.2, and excited in the TEoi5 mode, 
the transmission characteristics of the line are modified by 
the magnetic effect. The electric effect caused by the high 
dielectric constant DR can be assumed to be very small, 
since most of the transmission line is not only concentrated 
under the strip conductor but also almost orthogonal to the 
electric field of the DR [31] . The magnetic effect is an 
interaction between the magnetic field of the resonator and 
the magnetic field owing to the current in the strip 
conductor .This interaction can be regarded as mutual 
inductance Lm, and an emf is considered to be caused in the 
strip conductor. 

For this coupling the lateral distance between the strip 
conductor and resonator primarily determines the degree of 
coupling. Termination of the strip conductor and proper 
shielding also plays an important role in determining amount 
of coupling. 

Analysis 

Since the analysis of the configurations shown in the Fig. 
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5.2 is a formidable problem, an approximate method given by 
[32] and Komatsu [31] is used for analysis. 

The analysis is made on the basis of the following 
assumptions [31] . 

1 . The width of the strip conductor is much smaller than the 
diameter of the resonator. 

2. The transmission line carries only the TEM mode. 

3 . The DR can be represented by a conducing loop having in 
series an inductance Lr, a capacitance Cr and a resistance 
Rr. The resonance occurs at the frequency fo = iM (LrCr) . 

4 . Only the TEoi5 mode is strongly excited and field 
distortion caused by other modes is very small. 

5. The DR couples with transmission line through 

distributed mutual inductance Iji. 

With these assumptions, the structures of Fig. 5.2 can 
be represented by a resonance circuit magnetically coupled 
with a transmission line as shown in Fig. 5.3, Where Li, Ci 
and Ri are, respectively, equivalent circuit parameter for 
the transmission line terminated in Zi and Zo at load and 
feeding ends of it, respectively. Let Ir and Ii be the 
current flowing in the resonator and transmission line 
respectively. Applying KVL to this low frequency model we 
get 


V = IlRl + jwLlIl + jwLmlr + 


II 

IwCT 


0 


IrRr + jwLrlr + jwLmll + 


Ir 

jwCr 


eliminating Ir we get 


V=RiIi+jtt)LiIi+ 


(j^LbIi , Ii 

Rr + juLr + l/jwCr jwCl 


(5.1) 

(5.2) 


(5.3) 


The transfer impedance Zt, in the coupling plane is obtained 
as 


Zt 


y_ 

Ii 


. (i) 

+ R. 4. jmLr + l/Ji 3S •" 


3(3cr 


(5.4) 
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For perfect strip conductor Ri is very small and near 
resonance jwLi + 1/ (jwCi) may be neglected in comparison to 
terms containing Hence, the expression for Zt reduces to 

Zt = (5.5) 

Rr + - l/(LrCr)} 

Near resonance at u = wo + 5u transfer impedance may be 
written as 


Zt 


{(jJo+SO)) £iin 


Rr -1“ 




(5.6) 


Further simplification of Zt is obtained by neglecting 
(Sw)^ and (5(j)/wo terms 


Zt = 


2 2 

(Uo+SCj) Lm 


^ .{tjf+{5(j)|2Wo5a) -l/{LrCr) } 

+ Ci)o +oG) ^ 


(OoLm 


Rr . 2SW Lr 
+ 3 


0)o 


0)o 


Putting Rr = woLr/Qo in above equation we get 


Ci)oljini 

Lr . 25w Qo 
Qo 0)0 


Zt= Wo^^-Qo 


1 + 3X 


(5.7) 


Where 


X = 


2So) Qo 


(Jo 


Equation (5.7) indicates that the circuit shown in Fig. 
5 3 can be represented by a simple parallel tuned circuit as 
shown in Fig. 5.4, where L, C and E are obtained as 


40 



R=CJoQo- 


ll 


hr 




XjmCJo 


At resonance w = wo, so 


Zt (CJo) 


R = cjoQo- 


Lr 


(5.8) 


(5.9) 


here Qo is the unloaded quality factor of DR, 

For this parallel tuned circuit an important figure of merit 
is external quality factor Qe, which is the Q factor of the 
circuit when loading of circuit is only due to external 
Rext. Quantitatively it is given by 

Qe = Qa.-^^ (5.10) 


For our circuit, Rext symbolizes the internal impedance (Zo) 
of the source feeding to the line, and the load impedance 
(Zi) of the line. So, Rext = (Zi + Zo) . 


Qe = 


(Zi + Zo),, 
— R 


(Zi + Zo) 
OJo 


(Lr/li) 


(5.11) 


For known values of Zi, Zo and oo, the external quality 
factor Qe can be evaluated if we calculate (Lr/Lm) . The crux 
of the problem lies in finding the quantity (Lr/Lm) in terms 
of distance between DR and strip conductor. 

The DR acting in TEoiS mode can be assimilated to a 
magnetic dipole. Let Ir be the current flowing in the loop. 
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then, the voltage induced in the micros trip line due this 
current is given by 


e = -jWoLmIr (5.12) 

This induced voltage is also given by (Faraday's Law) 

e = -jwo/ioJgH.djs (5.13) 

where H is the magnetic field of the resonator and s, i.e. 
the surface of integration, is the surface perpendicular to 
the conductors passing into the substrate through which 
magnetic lines of H field crosses the strip conductor, 
from eqs.(5.12) and (5.13) we get 

L« = H.ds (5.14) 

Xr s 

The magnetic energy stored due to DR is given by 

Wm = (LrI?)/2 (5.15) 


From eqs. (5.14) and (5.15) we can obtain 


Lr _ 2 Wm 

Lm^ 4o{j'^H.ds}^ 


(5.16) 


Substituting the value of the expression from eqs. (5.16) 
into (5.11) we get 

Qe = (5.17) 

Oioiio {X^H.dfi}^ 


Now we see that for two configurations of Fig 5.2, only the 
surface integral will change with separation between the DR 
and the strip conductor. And also it will have different 
expressions for two configurations. Evaluation of XH.ds^ for 
both the configurations is given in appendix E, and the 
magnetic energy Wm was derived in chapter 3 . 
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5.3 Numerical Results 


5.3.1 CQtnparigQn with Experimental Data 

For the validation of present approximate analytical method 
presented above, computed Qe was compared with the 
experimental data [32] and shown in Fig. 5.5. Here, a 50-n 
microstrip line as a special case of suspended stripline, 
terminated its characteristic impedance at both ends is used 
for calculating Qe. Good agreement between calculated and 
experimental results justifies the use of the approximate 
analytical method for design purpose. 

5.3.2 Coupl ing between Suspended stripline and DR 

Fig 5.6(a) and Fig5.7(a) are showing the variation of Qe 
as a function of separation between the suspended stripline 
and DR, operating in S and Ka bands of microwave 
frequencies, respectively. There are three curves in each of 
the figures, for three different thicknesses of the strip 
conductor to show the effect of thickness on Qe, including 
the zero thickness case. 

5.3.3 Coupling between Inverted Microstripline and DR 

Fig 5.6(b) and Fig 5.7(b) are showing the variation of 
Qe as a function of separation between the inverted 
microstripline and DR, operating in S and Ka bands of 
microwave frequencies, respectively. There are three curves 
in each of the figures, for three different thicknesses of 
the strip conductor including the zero thickness case. 
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Chapter 6 


CONCLUSION AND FUTURE DEVELOPMENTS 

The dielectric resonator is studied in suspended substrate 
environment. The design parameters like resonant frequency, 
field configuration, energy distribution, coupling between 
two resonators and coupling between dielectric resonator 
with planner transmission lines are derived, in suspended 
substrate environment for millimeter wave applications. The 
main focus is on most commonly used TEoi5 mode of 
cylindrical resonator. The resonant frequency is calculated 
using an improved effective dielectric constant technique 
with well known dielectric waveguide model (DWM) method. 
Then the fields and the the energy distribution in all the 
regions of the structure is obtained. Using a magnetic 
dipole model for DR, the magnetic coupling is studied 
between two identical cylindrical resonators, for three 
important configurations. Coupling between cylindrical 
resonator and two important transmission line configurations 
like suspended stripline and inverted microstripline with 
finite thickness conductor, is studied using, simple lumped 
element model for DR and the line. 

Software is developed using above analysis, which can be 
used for CAD of microwave circuits like filters, 
oscillators, mixers etc. The results are generated for 
microstripline configuration as a special case of suspended 
substrate configuration, and are compared with experimental 
data available in the literature. This comparison validates 
the present analysis. 

Future Developments 

To get the full application of Suspended substrate 
environment, the hybrid mode analysis of dielectric 
resonator is required. Practical realization of the 
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microwave circuits like filters, mixers, oscillators, etc., 
using DR in suspended substrate environment is required to 
verify the configurations like off layered DRs and inverted 
microstripline with DR. The Theoretical results are very 
interesting for these cases, but due to non availability of 
the experimental data in literature for these 
configurations, we could not verify it. 

For upper band in millimeter wave range, the size of DR 
become very small, this problem can be solved by using 
Whispering gallery mode resonance in cylindrical resonators 
which allows the use of oversized resonators. 
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Appendix A 


Value Of Field Constants Bis 


/3i . tanh(aaha) .tanOih/2) - aa 
/3i . tanh (aaha) + aa ,tan(pih/2) 

[cos(Pih/2) + Bi ■ sinOih/2) ] 

e"C^*d^sinh(aaha) 

[cosOih/2) -Bi.sin(|3ih/2) ] 
[cos 02 h/ 2 ) - B4.sin02h/2] 

ocn - /32 ■ tanh(o:shs2) .tan(p2d3) 
^2. tanh(ash.s2) + as .tan(^2d3) 

B3. [cos (p2d3) - B4.sin(P2d3) ] 

e °^^d^sinh (ashs 2 ) 


Be = B 2 .B 7 


B7 


Jo (krd) 
Ko (ka.d) 


Be = B3.B7 


B9 = B5.B7 


where di = ha + h/2, d 2 = hs 2 + d3, and d3 = hsi + h/2. 
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Appendix B 


Expressions for Aijs and f sameS 


Aij Expressions: 


h/2 

A16 = X [cos(/3iz) + Bi.sinOiz)] dz 
-h/a 


, h/2 

A16 = S [sinOiz) - Bi.cosOiz)] dz 
-h/2 


h/ 2 

A18 = S sinh{as (z+hs2-h/2) } [cosOiz) - Bi.sinOiz) ]dz 
-h/ 2 


-h/2 2 

>137 = X [cos 02 z) + B4.sin02z)] dz 

-d3 


/ -h/2 2 

A37 = X [sinOaz) - B4.cos02z)] dz 

-d3 


hs 1 +h/2 

>137d =X [cos{p 2 (z-h 8 i -h)} 
h/2 


+ B4 . sin{p 2 {z-hsi -h) }] x 
[cosOaz) - B4.sin{[32z) ] .dz 
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tsamcS Expressions: 


t.a»e (h»2, d2, a.) :n [ Siph < 2ashB2 ) _ - 2«shs2] 


^ 2 a 8 d 2 . 

e .4at 


tsame (ha, dl , tta) ~ ^ttaha] 

aaadi 


e .4aj 


\ [sinh (2ashs2) + 2ashs2] 
t same (hs2, d2, tts) = ~ 


2 asd 2 ^ 
e .4as 


t same (ha, dl, Ota) =- 


[Sillh(20taha) + 20taha] 

20tadl . 
e .40ta 
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Appendix C 


Integrals involving Bessel Functions 
Ijo = r. jf (krr) . dr 
Iko = r .K? (kar) .dr 

d 

Icj = r.J^(krr). dr 
lek = X^ r.Kf(kar).dr 

d 1 

Ijk=X^^X ^r. [r -I- (2d+S) .sin(^)] J i (krr) Ki (ka.r" 
0 0 

where r' = [r^+ (2d+S)^ + 2.r. (2d+S) .sin(^) ] 


)drd^ 

1/2 
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Appendix D 


Evaluation of voliane integral XJ'/g.jI 

The volume in above integration is the volume occupied by DR 
and its supporting substrate in which polarization current 
exists. For edge coupled case we are considering second 
medium as dielectric material with dielectric constant csz 
and for suspended substrate case, this become air i.e. es 2 = 

1. While for off layered dielectric resonators this medium 
is always air. 

D.l M.ge Coupled case 

Here the integral can be splited in to different regions, and 
after substituting the value of polarization current from chapter 4 
gives 

X E.*> dv = jwoGo{ (ed-i)X (ELEi)dv + {esi-i)X (E7.E3)dv 

V V ^ — 

+ (es2-l)X^ (E8.E4)dv } (D.l) 

where gj denotes the electric field due to ith resonator in 
its jth region. 

The scalar product (Ej.gi) may be obtained in terms of 
the X and y components of electric fields. For this purpose 
we are considering the geometry as shown in Fig. D.l ,it is 
the bottom view of our structure (Fig. 4.2) which remains 
same for all three configurations, for off layered cases 
only we have to take care of z-coordinate transformation. 

According this figure we find that 

Ej.Ei = Ej El (cosy.sin^ + sinr.cos^) {D.2) 


where 

cosr = (y+y*) /r' 

sinr = x/r' 

2 


2 


sin0 = y/r 
cos^ = x/r 



Appendix E 


Derivation of Integral /l.dfi 

To incorporate the effect of thickness of the strip 
conductor, we have considered the induced emf on both sides 
surfaces of the conductor, and then taken average of it to 
get final expression of IH.ds, which is representing the 
induced emf on the conductor. On both sides of the conductor 
we will break the surface of integration s into parts 
according to the region number. The trick is same for both 
the cases i.e. suspended strip line and inverted 
microstripline. We are representing the surfaces as stop, 
and sbottom for top side plane of the strip and that towards 
the ground conductor, respectively. 

For suspended stripline case, we can write 

stop = SS + S6 and Sbottom =S7 + S8 (E.l) 

and, for inverted microstripline, we can write 

Stop = SS + S6 + 87 and sbottom = S8 (E,2) 

After deciding the surface of integration, now we can 
consider the Integrand i.e, H.dg. As we know that for TEoiS 
mode we have Hr and Hz components of H, but Hz component is 
normal to the surface vector so its dot product with the 
surface will be zero. Now only we have to consider Hr 
component. As shown in the Fig E.l the Hr.sin^ component is 
again normal to the surface vector so by same region it also 
vanishes. Hence, we have finally the integrant w.r.t. the 
Fig. E.l, as follows, 

(H.ds)i = Hri.cos^ dx dz {E.3) 

where the subscript i denotes the region number in which the 
integration is to be performed. 

From the geometry of Fig. E.l we can write 

COS0 = b/{x% and r = {x^ (E.4) 


After sustituting the value of cos<f> and that of the r in to 



field expressions of H, we get the final expressions of 
integration in different regions. Then put these values into 
eqs . (E.l) and (E.2) , for suspended stripline and inverted 

microstripline respectively. 

Finally 

Jl.dfi = [| {X Ii.d£}stop| + I {X H.ds}sbotto,n|3/2 {E.5) 

For both the cases, 


{XH.ds,}gg= Be.exp (-aadi) sinh{aaha) .Ih/ka 

{Xii.ds}g^=B8 [cos (/32d3) -cos (|32h/2) 

-B4{sin (j32d3) -sin(/32h/2) }] Ih/ka 

But for suspended stripline: 


{E.6) 


{Xii-ds}_„= - B9. exp (-otsd2) sinh(ashs2) . Ih/ka 

o o 


{Xil.dfi} _= B 7 [cos Oih/2-t) -cos Oih/2) 

S' o 

- Bi{sin(/5ih/2) + sin{[3ih/2-t) }] Ih/lca 

(E.7) 


and, for inverted microstripline: 


{XS.d£}„^= - 2 BiB 7 sin(/3ih/2) Ih/ka 
S o 

{XH.ds,} = - B 9 . exp (-asda) sinh(ashs2-t) . Ih/ka {E. 8 ) 

In the above expressions, t is the thickness of the strip 
conductor. The coupling plane between DR and strip-conductor 
is at li distance from feeding end and at I2 distance from 
terminal load where the length of the line is (I1+I2) and 
width is w. This is shown in the Fig E.l. The integral Ih is as 
follows : 


Ih = X^^jKi{ka(x^+ b^)^^^}.b/{(x^+ b^)^^^}dx 


(E.9) 
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